From Shock Waves to Brownian Motion and 1/f- Noise in Gas 



O 

o 

(N 



I 

o 



> 
o 

00 

in 

in 
o 
i> 
o 



X 



Yuriy E. Kuzovle\EI 

Donetsk Institute for Physics and Technology (DonPTI NASU), 83114 Donetsk, Ukraine 

A formally exact relation is derived which connects thermodynamically non-equilibrium evolution 
of gas density distribution after its arbitrary strong spatially non-uniform perturbation and evo- 
lution of many-particle correlations between path of some marked particle and its surroundings in 
equilibrium gas. This relation directly confirms significance of the many-particle correlations even 
under the Boltzmann-Grad limit and thus validates the earlier suggested revision of kinetics. 
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I. INTRODUCTION 

The "law of large numbers", first discovered by 
Bernoulli approximately 300 years ago up to now is 
world outlook paradigm of physical applications of prob- 
abilities. Meanwhile, at the end of his work [l| Bernoulli 
himself emphasized that if equivalent observations were 
perpetually continued (and at last probability turned to 
certainty) then one would conclude that all things obey 
exact rules or even fate. In the field of probabilistic inter- 
pretation of quantum mechanics such conclusion means 
"nonlocality" not less impressive than the famous EPR 
paradox. Naturally, one may ask who is counting out- 
comes of our observations (e.g. castings a die) and makes 
so that their relative frequencies have definite limits equal 
to speculative "probabilities" (e.g. all equal to 1/6) ? 

The standard answer is that those godlike role belongs 
to "independency" of the outcomes, which claims that 
probability of their series is product of elemental prob- 
abilities. But this is tautology, not the answer, because 
"elemental probability" loses meaning in absence of the 
limit. A good explanation of this logical fact can be found 
in 0. Hence, presence (or absence) of the limit is noth- 
ing but postulate concerning probabilistic measure in a 
space of observation sequences, and no more. 

Unfortunately, in statistical physics built on micro- 
scopic dynamics the only a priori probabilities at our 
disposal are probabilities of its initial conditions. In 
such the theory, as Krylov emphasized almost 60 years 
ago 0, we should distinguish between actual dynamical 
(cause-and-effect) dependencies or correlations and sta- 
tistical ones. Moreover, sometimes they controvert one 
to another. The determinism of dynamics can produce 
infinitely long-living actual correlations and thereby ar- 
range definite limits for relative frequencies being per- 
ceived as fast decay of statistical correlations and "in- 
dependency" of observations (similarly, ideality of a die 
is incarnate actual correlation ensuring "independency" 
of its castings). Many such examples are known from 
mixing (time-reversible chaotic) nonlinear dynamics with 
finite number of degrees of freedom. 
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And, vice versa, deficiency of actual dynamical de- 
pendencies may result in arbitrary long statistical cor- 
relations. For example, if earlier collisions of a particle 
(or charge carrier) with gas of similar particles (or with 
phonon gas) have no infiuence on later ones, then such 
system is indifferent to what rate (relative frequency) of 
collisions does happen (therefore it would be vain to wait 
till "probability turns to certainty" ) . The result is scale- 
less fiuctuations, possessing 1/f -type spectrum, in col- 
lision rate of a test particle and thus in its diffusivity 
and mobility ( 1/f -noise) [1, d, 0] (for more details and 
examples see also [1, H, ll]). 

As it was shown in 01 (see also [1, i [I^l), 

in such 

case one has no rights to introduce a priori "probabilities 
of collisions" , all the more, postulate "molecular chaos" . 
Instead, we should honestly investigate dynamical evo- 
lution of an initial statistical ensemble, taking in mind 
Krylov's conjecture that under mixing dynamics, gener- 
ally, relative frequencies of some phenomenon along par- 
ticular phase trajectories in no way are connected with 
some its a priory probability @ . 

In case of gas, if wishing to investigate Brownian mo- 
tion of its particles, we have to solve the BBGKY equa- 
tions starting from spatially non-uniform statistical en- 
semble with certainly known (non-random) position of at 
least one (test) particle. An exact solution is wittingly 
impossible, but we can presume the gas not too dense 
and resort to approximative shortened description of col- 
lisions in terms of Boltzmannian "collision integrals". 
However, analysis of this approach in revealed that 
the Boltzmann's "molecular chaos" hypothesis is incom- 
patible with spatial non-uniformity of the ensemble. As 
the consequence, in place of the single Boltzmann equa- 
tion the theory births infinite chain of kinetic equations. 
They operate with usual one-particle distribution func- 
tion plus infinite set of specific many-particle distribution 
functions which represent clusters of (real or virtual) col- 
lisions or close encounters of particles and together de- 
scribe fluctuations in the rate of collisions 0, B Q ■ 

In view of innovatory character of these specific func- 
tions, it seems principal to make sure that they have 
adequate prototypes beyond approximations. Below we 
present a proof based on exact relations between such 
functions and a relaxation law for arbitrary strong per- 
turbations of gas density. 
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II. START STATISTICAL ENSEMBLE 

First, we will consider statistical ensemble of N par- 
ticles which are contained in a volume fl and possess 
short-range repulsion of one from another, under two 
successive formal limits: firstly the "thermodynamical" 
one, N ^ oo , Sl^oo, = N/il. ^ const, and sec- 
ondly the "Boltzmann-Grad limit", vq oo, cr — > , 



• const ( A is mean free path and a prop- 



erly defined cross-section of inter-particle collisions). 

Let initially, at time moment t = , our system is 
described by the A'^-particle distribution fmiction 

N 

p(-)(r) = C(*"V('='J)(r) exp[-^t/(r,)/T] , 

p(e9)(r) ^ c'^'^<}) exp[-iJ(r)/T] (1) 

Here F = {Fi,...,FAr} with Tj = {vj ,pj} being vari- 
ables of j-th particle, H denotes full A^-particle Hamil- 
tonian, p^^"^-* is corresponding canonical equilibrium dis- 
tribution , V{r) is a bounded function, and constants 
C^'^'?) and C(™) ensure the normalization 



j =1 j =1 



(2) 



We will be interested in consequent evolution of the 
system, that is in distribution function 

exp(tL)p(-)(F) , 



L 



N 

E 



[dH/drj)d/dpj - (dH/dpj) d/dr^ 



(3) 



where L is the Liouville operator corresponding to 
Hamiltonian H . Clearly, p^™^ and p{t) represent sta- 
tistical ensemble which is non-equilibrium because of 
non-uniform perturbation of particles density 



(4) 



Alternatively, one can think that p{t) comes from equi- 
librium in presence of an external potential V{r) which 
however is abruptly switched off at t = Q . 

It is convenient to characterize density perturbations 
by function 

4){r) EE exp[-y(r)/r] - 1 
Then, introducing functional 

N 



^{0} ^ / p'^^'^\T) n[l + 0(r,)] dF 



one can easy verify the equalities 



1.(0, r) = [l + <^(r)] 



(5) 

(6) 
(7) 



In order to perform the thermodynamical limit, we will 
assume that as a whole the perturbation is finite: 



\<j){r)\ dr < oo 



Then the limit of the functional ([5]) looks as 



(8) 



k. r 

= 1 + ^!^ / Xl<i>{r,)dr, = 

k=i ' j=i 

OO P k 

J„ 1 'J A 1 



where Fk{...) = Fk{ri, ... ,rk) are usual non-normalized 
(or, in other words, "normalized to volume") /c-particle 
spatial distribution functions (or, in other words, "corre- 
lation functions") of strictly equilibrium gas. In particu- 
lar, Fi(ri) = 1. 

Functions Tfe(...) = Tfc(ri, r^) formally relate 
to Fk like cumulants, or semi-invariants, of a random 
field relate to its statistical moments: T2(ri,r2) = 
F2{ri,r2) - Fi{ri)Fi{r2) = F2{ri,r2) - 1 , and so on. 

In fact, of course, all F^ and obey translation 
invariance, i.e. depend on inter-particle distances —rj 
only. At that, any of the cumulants T^; — > when some 
of its arguments — rj goes to infinity as measured in 
units of characteristic radius rg of particles interaction 
(taking in mind 3-D gas, one can define rg e.g. by means 
of cr = tttq ). Therefore, if the function (j>{r) is smooth 
enough, on scales of order of tq , then 



fe 



rj)drj 



Ik To 



3(fe-l) 



0'=(ri)dri (10) 



at /c > 1 with 7fc being "virial coefficients" . Under the 
Boltzmann-Grad limit, rg — > and the "gas parameter" , 
defined e.g. hy g = {AirrQ /?>) vq , also tends to zero, 
(7 — > . Consequently, gas becomes ideal from the point 
of view of thermodynamics: 



T{(t)} — > exp 



vo I (j){r)dr 



(11) 



^0 



l + <t>{r) ^ cxp[-l/(r)/T] (12) 



III. 



RELAXATION, TIME SYMMETRY AND 
EQUILIBRIUM-NON-EQUILIBRIUM 
RELATIONS 



However, from the point of view of kinetics our gas 
remains quite non-ideal. Let us consider relations which 
substitute ^ and (O at t > and connect, on one 
hand, non-equilibrium process of relaxation of relative 
gas density perturbation, v{t,r)/vQ — 1, and, on the 
other hand, Brownian motion of particles in unperturbed 
equilibrium gas. 
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Taking into account that the Liouville operator L 
commutes with any equiUbrium distribution, L p^'^'^^ f = 
p(<^i)ijf (where / is "arbitrary function" of Tj), and 
combining formulas ((T|), ^ and we can transform 
the latter to 



1 r ^ 



(13) 



N 

E 



N 



<5(r-r,)exp(ti) ll[l+^{r,)] 
j=i 



Let us make the change of integration variables here: 
Pj " Pj ■ Since the Hamiltonian is even function 
(quadratic form) of the momentums while L is odd 
with respect to them, the only result of this operation 
is change of sign of time argument before L : 



1 



N 

E 



N 

n p 

N 



(eg) 



(T) X 



(14) 



Sir-r,)cxpi-tL) + cj,{r,)] 



which clearly demonstrates time symmetry of gas evolu- 
tion: i^{—t,r) = i'{t,r). 

Next, recall that operator exp ( —tL) has twofold 
sense: if acting onto distribution functions it describes 
their time-reversed evolution, but if acting onto system's 
variables (or any function of them) it describes their evo- 
lution in real direct time. In particular, 



exp(-iL) ^{r,) = 0(rj(i,r)) 



(15) 



where expression rj{t,T) means (vector of) coordinates 
of j-th particle at time t considered as a function of ini- 
tial state r of the whole system taken at t = (i.e. the 
set rj {t, r) is solution to system's Hamiltonian equations 
starting from rj{0,r) = rj , pj{0,T) = pj) . 

Hence, obviously, we can rewrite p4p in the form 



p('=«)(t,r|r) 



N 



{r,)]dT, (16) 



where /9^'^''^(t, F | r) is conditional distribution function 
of system's state F at current time t corresponding to 
condition that initially (at time zero) one of particles was 
positioned exactly at a given point r . Formally, 



{t,T\r) - exp(tL)p('='')(0,F|r) , 



(0,F|r) =J7p('=^)(F)^ 



N 



N ^ 



5{r-r,) , (17) 



where factor fi/iV ensures normalization to unit: 
/ p'^'^''^ (t, F I r) dF = 1 , and we already keep in mind the 
coming thermodynamical limit. 



To complete the limit, notice that all of A'' terms of 
(fT7)) give equal contributions to integral in In anal- 
ogy with transition from ([5|) to (l9|), the result is 



1 



dr' [l + (t){r')]{Wiit,r'\r) + 



oo ^ „ k 

+ E9"/ F>^+i(t^^Vi^ -^^klr) Y[<l>{r,)dr,} (18) 
fe=i ■ j=i 

Here Wi{t,r'\r) means probability density of finding 
at point r' at time t the same test particle what started 
from the point r . This function is normalized to unit: 



/ Wi{t,r'\r)dr' = 1 



(19) 



In other words, this is just the distribution of Brownian 
displacement of test particle which was under investiga- 
tion in [§] (and much earlier in 0, Si)- 

Functions Fk+i{t,r'\ri, ... ,rk\r) describe joint prob- 
abilities of finding the test particle (started from r ) at 
point r' and simultaneously some k other particles at 
points Tj . These functions are symmetric with respect 
to (j = l...k) and definitely normalized with respect 
to r' as well as r : 

J Fk+i{t,r'\ri,...,rk\r)dr' = Fk{t,ri, ... ,rk\r) (20) 

J Fk+i{t,r'\ru-,rk\r)dr = F^+i (n, ... , r^, r') (21) 

On right-hand side in (PH)) we see non-normalized (i.e. 
"normalized to whole volume") distribution functions 
(correlation functions) characterizing probability of lo- 
cation of k particles at points rj under condition 
that one more particle had started from r . Functions 
Ffc+i(ri, ... ,rk,r') on right-hand side of (pij) are usual 
static (fc + l)-particle distribution functions mentioned 
in previous Section. 

Of course, after the limit all these functions acquire 
translation invariance: Wi{t,r'\r) = Wi{t,r' — r) , 
etc. With this notion it becomes rather obvious that 
identities (PT|) guarantee "conservation of particles" : 
/ v[t,r)dr = J i/{0,r)dr. 

Initial conditions for these functions also are clear: 

Wi{0,r'\r) = 6{r' -r) , 
Fk+i{0,r'\ri, ... ,rk\r) ^ S{r' - r) Fk+i{ri, ... ,rk,r) , 
Fk{0,ri,... ,rk\r) = F^+i (ri, ... , r^, r) 

(22) 

Correspondingly, at i = ^ reduces to ^ with 

In principle, (|18p can be qualified as example of "gen- 
eralized Onsager relations" or "generalized fluctuation- 
dissipation relations" which follow from time reversibility 
of microscopic dynamics and establish connections be- 
tween dissipative nonlinearity and statistics of transport 
processes, in particular, statistics of equilibrium noise 
(see e.g. [lllli^l and references therein). 
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IV. DYNAMICAL MULTI-PARTICLE 
CORRELATIONS AND THEIR CUMULANTS 

The distribution functions what appear in (|18|) mix up 
i) average values, ii) static correlations between positions 
of particles due to their potential interaction, and be- 
sides iii) specific dynamical correlations between current 
relative disposition of mutually colliding (or merely en- 
countering) particles and their previous displacements . 
Next, we would like to separate those contributions and 



especially consider the third of them. 

In order to extract joint dynamical correlations in 
a pure form, we should introduce cumulants, or semi- 
invariants, similar to cumulant functions Tfc(ri, r^) 
in With this purpose, first, notice that expression 
inside braces in (fT8|) has the same formal structure as 
expression ( x exp (<j)y)) , where x and y are two ran- 
dom quantities, cf) is dummy parameter (test parameter), 
and angle brackets denote averaging procedure. Second, 
apply the well known general formula (see e.g. ilSj) 



( a; exp ((/) y) ) = ( exp (0 y) ) 



{x) +(j){x,y) + — {x, y,y) + Y^ {x, y, y, y) + 



(23) 



with {x,y, ... ,y) being joint cumulants in Malakhov's notations [13| . 
In our context, similar transformation as performed in (|18p yields 



C5C ^ „ ^ 

W,{t,r'\r) + J2'^ / Fk+i{t,r'\n,...,rk\r)l[<l>{r,)dr, 
k=i - ■' j=i 

^ + Fk{t.r,,...,rk\r)l[c^ir,)dr, 

k=i ' j=i 

CXD ^ „ k 

Wi{t,r'\r) + Y.T^ / Wk+i{t,r'\ri,...,rk\r)\{<jy{r,)drj 



(24) 



J^{t,0|r} X Y\J{ty\4,\r} 



Here new functions Wfc+i (<, r' |ri, ... , r^jr) , in place of Ffe+i (t, r' |ri, ... , rfe|r) , are just the cumulant functions 
(analogous to (x,?/, ...,?/) in (|23p ) what correspond to refined dynamical (two-time) inter-particle correlations. 

Further, notice that because of the normalization relations ([T9|) and ([20| result of integration of the above expression 
over r' coincides with functional 0|r} defined by the upper square brackets in Therefore, firstly. 



y" yV{ty\<\)\r} dr' ^ \ , thatis ^ VFfc+i(t, r 'In, ... , rfc|r) dr' = , (25) 

with functional yV{i, r '|(/) |r} defined by the lower square brackets in (|24p . Secondly, formula (HH]) takes the form 

v{t,r) T{t,(^\r} 



1 



(r')>V{t,r'|</>|r} dr' 



Notice also that, because of the normalization relation (|2ip . integration of ([24]) over r results in identity 

1 (5J^{(/.} 



y" T{t,^\r} >V{i,rV|r} dr = 



1^0 5(t>{r') 

Besides, combining (j24[) with initial conditions (|22p we obtain one more property of cumulants Wk+i '■ 

W{0,r'\(f>\r} ^ 6{r' -r) , thatis Wk+i{0,r'\ri, ... ,rk\r) = , 

thus confirming that indeed they represent correlations of very dynamical nature. 

Functions Fk{t,ri, ... ,rk\r) also can be replaced by corresponding cumulants if we write 



T{t,4)\r} = exp ^ / Tfc(t,ri,...,rfc|r) Jl 0(rj)drj , 
fc=i ' j=i 



(26) 



(27) 



(28) 



(29) 



with Tfc(t, ri, ... , rfe|r) being direct analogue of above introduced static cumulant functions Tfc(ri, r^) 



5 



Let us discuss and compare two just now introduced 
sorts of cumulants. Last of them, Tk{t,ri, ... ,rfc|r) , at 
k > 1 describe conditional joint correlations between k 
particles under condition that some one more {k + l)-th 
(test) particle earlier, at t = , was localized at point 
r . Certainly, such poor condition can not change short- 
range character of the correlations, and the latter must 
decay at same inter-particle distances jr^ — Tj| > tq as 
unconditioned correlations. Then, quite as in (jlOp . 



J Tkit,ri,... ,rk\r)f[(j){rj)dr 



^o'""'^ / 7fc(i,n |r )</)'= (ri)dri (30) 



All the more, such condition can not influence upon 
mean gas density. Therefore it is natural to state, in 
agreement with ((22|) . that difference between the condi- 
tional one-particle distribution function, Ti(i,ri|r) = 
Fi(t, ri\r) (expressing conditional gas density in units of 
vq), and the unconditioned one, Ti(ri) = Fi{ri) = 1 , 
at t > is not significantly greater than at t — . We 
can write this in the form 



y"[Ti(i,ri|r)-Ti(ri)](/.(ri)dri 



< 



(31) 



< 



[F2(ri,r)-l]0(ri)dri 



where (j) is characteristic magnitude of 4'{r) . 

Formal substantiation of these statements in princi- 
ple could be obtained from the Bogolyubov-Born-Green- 
Kirkwood-Yvon (BBGKY) equations 



dt 



n „ 

LWi^(^')+z.oE hn+iFif,dT,,+^ (32) 



for the full phase space n-particle distributions F, 
= Fi'^''^(0,ri|r2, ... ,r„_i|r) connected to our spatial 
distributions by means of 

W,{t,r^\r) = j Ft''\t,r,\r)dp, , 

Ffc+i(i,ri|r2, ... ,rfe+i|r) = (33) 
^ fc+i 

- / F,(^\Ht,ri|r2,...,r,+i|r) Hdp, 

along with ([^0]) . and subjected to initial conditions 

Fi'^«)(o,ri|r2,...,r„^i|r) = 

n 

= (5(ri-r)F„(ri,...,r„) n^o(p,) , (34) 

with "number one" assigned to the test particle, L*^"^ 
standing for n-particle Liouville operator, 

L,, = VU{r, - r,) • {d/dp, - d/dpj ) , 



U{R.) being interaction potential, and Go{p) equilib- 
rium Maxwellian momentum distribution. Clearly, this 
is a task about Brownian motion of the test particle in 
equilibrium gas. But, unfortunately, such BBGKY equa- 
tions still remain almost uninvestigated beyond rough- 
ening like the Boltzmannian kinetics or its improvement 
suggested in 0, § (otherwise essays like the present one 
would be unneeded). 

Other conditional cumulants, Wk+i{t, r' |ri, ... , r^jr) , 
seem much more interesting objects because describe 
joint {k + l)-order correlations between k "ordinary" 
particles and besides the test particle itself. The latter 
is peculiar since it is the only particle we know where- 
from it came to its current position. Therefore, in fact 
Wk+i represent joint correlations between its preceding 
displacement r' — r { "Brownian path" ) and its current 
surroundings. Moreover, such extra correlations what 
do not reduce to mere gas density fluctuations already 
taken into account by Tj, . Thus we must conclude that 
cumulants Wk+i{t, r' \ri, ... ,rk\r) represent fluctuations 
in the line of collisions of the test particle, and in this 
sense, in contrary to , they are specific dynamical 
correlations. Correspondingly, coordinates rj belong to 
those neighbors of the test particle which currently are 
or recently were or soon will be involved into collisions 
with it (or even those which would be in collision if some 
other neighbor did not interfere with it). 

As joint cumulants in general, Wk+i{t,r' |ri, ... ,rk\r) 
vanish if one (or more) of its constituents become sta- 
tistically independent on others. The independency, in 
its turn, takes place when any of the relative distances 
\rj — r'\ exceeds some characteristic spatial scale • 
However, the above reasoning prompts that this scale is 
not the interaction radius tq , figuring in pH)) and (PU)) . 
but "radius of collision" which is sooner the mean free 
path A than rg . 

Indeed, clearly, information about a collision propa- 
gates to as long distance apart it as A , if not greater. 
Nevertheless, A would be overestimate of Iq , since that 
is propagation along the direction of relative momentum 
Pj ~p' (and velocity) only. At two perpendicular dimen- 
sions the collision occupies merely a region w ct ~ Tq . 
Hence, total "volume of collision" (and volume covered 
by related correlations) is ~ <tA = I/vq. Correspond- 
ingly, it is reasonable to estimate mean (averaged over 
all directions) "radius of collision" (and radius of related 

1/3 

correlations) Iq as Iq = l/'^o • -'^^ other words, Iq 
is nothing but typical distance between neighboring gas 
particles. 

For preliminary characterization of cumulants M^fc+i 
it is useful to notice also that both the initial conditions 
(P5| and equalities ([^5)) are manifestations of the men- 
tioned general property of cumulants: since at t — po- 
sition of the test particle is non-random (certainly equals 
to r ), any joint cumulant it enters must disappear (sim- 
ilarly, all ( cc, y, ... , y ) = if a; is a constant), while (1^5)) 
as well reflects non-randomness of event "test particle is 
somewhere" . 
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V. THE BOLTZMANN-GRAD GAS 



Next, consider cumulants Wk+i under the Boltzmann-Grad limit, when their role becomes especially distinct. 
At g ^ , the estimates pU)) . (1501) and, very importantly, (pij) . lead to the asymptotic 



!F{t, (j) \r} — > exp 



t^o j (j){r)dr 



As the consequence, firstly, identity (P7|) simplifies to 

j W{t,r'\(j)\r} dr = I , that is J Wk+i{t,r'\ri, ... ,rk\r) dr = 
Secondly, our main formula ^TE\i . or turns into 

^^^^-1= /(/)(r')>V{t,r'|0|r} dr' 
Thirdly, definition (|24p of the cumulants Wk+i{t,r'\ri, ... ,rk\r) degenerates to 

Wi(i,r'-r) + ^ i f Fk+i{t,r'\ri,...,rk\r)Y[^{r,)dr, = 

fc =1 ■ J =1 

= I Wi{t,r' -r) + f2^^ J Wk+iit,r'\ri,...,rk\r)Y[^{rj)dr, i exp 



93(r) 



where we introduced new dummy function (series expansion variable) (p{r) = (f'i'f) ■ 

Fourthly, functional W{t, r '|0 |r} should be reconsidered in the course of the Boltzmann-Grad limit, 

oc ^. „ k 

W{t,r'\cb\r} = Wi{t,r' -r) + lim V ^ / r 'In, ... , rfc|r) ff 0(r,) dr, , 



(35) 



(36) 



(37) 



(38) 



(39) 



as a result of competition between infinitely growing fac- 
tor j/q , on left hand, and infinitely shrinking radius Iq 
of decay of Wk+i{t,r'\ri, ... ,rk\r)'s via. \rj~r'\ de- 
pendence, on right hand. 

Among three qualitatively different variants of the 
competition, i.e. vqIq — > or vqIq —> const or 
^0^0 ~^ oo J in fact the middle variant is acceptable only, 
since, from the point of view of the density relaxation, 
vit, r)/i>o , the third variant would mean that none ini- 
tial density perturbation can be considered as a weak one 
(even if it is infinitesimally small) while the first variant 
would lead to purely linear relaxation of the perturba- 



tion (even if it is arbitrary large). Hence, once again we 

1 /3 

come to conclusion that Iq ^ l/^^o (and without loss 

1 /3 

of generality we can put on = l/'^o )■ 

In order to properly perform the limits in p9p . one 
should take in the mind identities (and then ([M]) as 
well) which are very important since say that results of 
integrations in (|39p must represent full derivatives with 
respect to r' (and finally also with respect to r). 

Then, under assumption that (pir) is not "too arbi- 
trary" but sufhciently smooth field, quite natural expec- 
tation for the limit of (fc -t- l)-th term of (|39p reads 



Wk+i{t,r'\ri,...,rk\r)l[q)ir,)dr, = - . I,+,{t,r' - r) cj,'^ {r') (40) 
j =1 

with vector functions Ik+i{t,R) formally just defined by this expression. Substituting it into and then to 
(157)1 and twice applying integration by parts, we obtain 

^^-1- fw,it,r'~r)<j>{r')dr' + f]-l— fwk+iit,r'~r)<f>>'+^{r')dr' (41) 



with one more set of (scalar) functions defined by 



Wk+i{t,R) = -dIk+i{t,R)/dR 
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(42) 



According to these definitions and again to identities (I^S)) 
and (p8| . with addition of ([36)) . the new functions must 
obey restrictions 

J Wk+i{t,R) dR^ , Wk+i{0,R) = , (43) 

J Ik+i{t,R) dR = , Jfc+i(0,i^)=0 (44) 

We passed over in silence that in fact transition from 
left-hand to right-hand side of needs in one more 
property of cumulants Wk+i{t, r'\ri, ... ,rk\r) , namely, 
that they remain bounded under the Boltzmann-Grad 
limit. In terms of corresponding distribution functions, 
Fk+i{t,r'\ri, ... ,rk\r) < const . Orally, magnitude of dy- 
namical inter-particle correlations does not increase when 
their mean radius Iq decreases down to zero. From the 
point of view of the BBGKY equations, hardly this state- 
ment gives rise to doubts. What is for a scale of smooth- 
ness of 0(r) needed in (|40| when Iq , possibly, it 
must stay at the level of mean free path A . 

VI. DISCUSSION AND RESUME 

It seems quite obvious that just now considered distri- 
bution functions Fk+i{t,r'\ri, ... ,rk\r) , corresponding 
cumulant functions Wk+i{t,r'\ri, ... ,rk\r) and their 
children w k+iit, R) are in very close conceptual and 
mathematical connections with the special full phase 
space distribution functions An{t, R,Vi, ... ,Vn) , consid- 
ered in (tJ in the framework of the "coUisional approxima- 
tion" to BBGKY equations (in Q they were renamed to 
Fn{t, R, Vi, ... , Vn) ), and their children like Wn{t, R) in- 
troduced and investigated in [9^]. At that, equalities ([^5]) 
and (UD) show that, in opposite to Wi{t,R), functions 
Wk+i{t, R) can not be treated as spatial probability dis- 
tributions. The lowercase letter just accentuates this dif- 
ference between w k+i{t, R) and true (non-negative) dis- 
tribution functions Wn{t,R) . By implication, however, 
Wk+i{t,R) and Wk+i{t, R) reflect similar properties of 
statistical ensemble. More careful analysis, with the help 
of the BBGKY equations, shows that sooner the first 
represent time derivative of the second. 

Detailing of these connections is interesting but non- 
trivial subject for separate work. At present, the thing 
what is most important for us is that in essence both 
the old and new of the enumerated mathematical objects 
represent the same dynamical multi-particle correlations. 
This notion, as combined with formally exact relation 
(PS)) (or its limit form dH])), leads to principal conclu- 
sion that the specific multi-particle distribution functions 
An{t, R,Vi, ... ,Vn) , half-intuitively introduced and ar- 
gued in as autonomous complementary (in addition to 
the one-particle distribution) characteristics of spatially 



non-uniform statistical ensemble, now appear as formally 
well-grounded concepts of statistical kinetics. 

The salt of such statistical characteristics of gas as 
An{t, R,Vi, ... ,Vn) or Wk+i{t,R) is that they were ir- 
retrievably lost by conventional kinetics based on the 
Boltzmann's "Stosshalansatz" or, in other words, "molec- 
ular chaos" hypothesis. Due to it, in the Boltzmannian 
kinetics (at least under the Boltzmann-Grad limit) the 
one-particle distribution function ascended the throne as 
the only and exhaustive characteristics of gas state, since 
none inter-particle correlations have survived. But now, 
at last, we have visual and, in my opinion, rather con- 
vincing demonstration of falseness of such kind of kinetic 
theories. 

Concretely, to resume, we demonstrated existence and 
significance of dynamical multi-particle correlations what 
follow from the above derived chain of relations (HH]) — > 
(I26l)^(l37l)+(l39l)^(l4l|). These relations connect, on 
one (left) hand, 

coefficient functions in series expansion of space-time- 
varying density of spatially non-uniform and thus ther- 
modynamically non-equilibrium gas over its initial den- 
sity perturbation, and, on the other (right) hand, 

cumulant correlation functions which describe creation 
and evolution of joint many-particle correlations in ther- 
modynamically equilibrium gas which is spatially uni- 
form in thermodynamical sense but at the same time 
non-uniform in statistical sense owing to "informational 
perturbation" of gas by information about initial position 
of one of its particles (termed "test particle" ) . 

On the left, initial density perturbation can be arbi- 
trary large, therefore generally it causes strongly non- 
linear relaxation process, may be even producing shock 
waves. This means that on the right-hand side all many- 
particle correlations really exist (remaining essentially 
different from zero even under Boltzmann-Grad limit). 

At it was shown earlier, first in 7] (see also [8]) and 
then in [lo| , these correlations are not things in them- 
selves but in fact govern a true evolution of the one- 
particle distribution function (referred to the "test parti- 
cle" ) and thus determine statistics of self-diffusion (Brow- 
nian motion) of gas particles, at that involving possibility 
of what can be characterized as "violation of the law of 
large numbers" and "1/f fluctuations in particles diffu- 
sivity and mobility" . Immediate cause of the loss of these 
correlations at very beginning of kinetics was the way of 
thinking when one identifies actual independency of par- 
ticles or events, in the sense of dynamics, and their formal 
independency, in the sense of mathematical probability 
theory (see Introduction). The "molecular chaos" hy- 
pothesis is natural product of such thinking about "colli- 
sions" of particles. However, this hypothesis immediately 
becomes inadmissible if "collisions" are presumed to sat- 
isfy such principal property of dynamics as conservation 
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of phase volume and consequently conservation of prob- 
abilities (for more explanations see 0, [1] and references 
therein) . The pay for so accurate definition of collisions 
is absence of definite limits for their relative frequencies, 
i.e. absence of a priory "probability of (one or another 
sort of) collisions" at all (paraphrasing Bernoulli's words 
[3], probability never turns to certainty). That is why 
"molecular chaos" hard stays unproved (see article "On a 
derivation of the Boltzmann equation" by O.E. Lanford 
in book and related comments in [8|), and instead 
of the single Boltzmann equation we have to deal with 
infinite chain of equations for specific distribution func- 
tions which together represent time-scaleless fluctuations 
in "probability of collisions" . 



Up to now, such revision of kinetics was concerning 
statistics of self-diffusion in equilibrium gas only 0, H, 0] . 
But if continued to non-equilibrium situations it must 
somehow touch also kinetics in its usual sense of the fluid 
mechanics. Intuition would like to suggest that probable 
changes in this field must be not revolutionary. If that 
is the case, then it is interesting, for instance, to clar- 
ify what information about the Brownian motion follows 
from the mentioned relations if one substitutes there re- 
sults of conventional kinetics, e.g. those coming from the 
Boltzmann equation (though, to the best of my knowl- 
edge jl4l |. none exact spatially non- uniform solution to 
this equation is known) . In any case, we have found new 
incentives to critical review of kinetics. 
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